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Abstract 



A method based on order completion for solving general equations is 
presented. In particular, this method can be used for solving large 
classes of nonlinear systems of PDEs, with possibly associated initial 
and/or boundary value problems. 



"... provided also if need be that the notion of a solution 
shall be suitably extended ..." 

cited from Hilbert's 20th Problem 



1. Preliminaries 

Recently in [3], systems of nonlinear PDEs composed of equations of 
the general form 

(1.1) F(x,U(x), . . . ,DPU(x), . . . ) = f(x), i60 C W 1 
were solved on domains Q that can be any open, not necessarily 
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bounded subsets of IR n , while p e N n , |p| < m, with the orders m G N 
of the PDEs arbitrary given. 

The unprecedented generality of these nonlinear systems of PDEs comes, 
above all, from the class of functions F which define the left hand 
terms, and which are only assumed to be jointly continuous in all of 
their arguments. The right hand terms / are also required to be con- 
tinuous only. 

However, with minimal modifications of the method, both F and / 
can have certain discontinuities as well, [3]. 

Regardless of the above generality of the nonlinear systems of PDEs 
considered, and of possibly associated initial and/or boundary value 
problems, one can always find for them solutions U defined on the 
whole of the respective domains Q. These solutions U have the blan- 
ket, type independent, or universal regularity property that they can 
be assimilated with Hausdorff continuous functions, [1,4-6]. 
It follows in this way that, when solving systems of nonlinear PDEs of 
the generality of those in (1.1), one can dispense with the various cus- 
tomary spaces of distributions, hyperfunctions, generalized functions, 
Sobolev spaces, and so on. Instead one can stay within the realms of 
usual functions. Also, when proving the existence and the mentioned 
type of regularity of such solutions one can dispense with methods of 
Functional Analysis. However, functional analytic methods can pos- 
sibly be used in order to obtain further regularity or other desirable 
properties of such solutions. 

The mentioned generality of the equations solved and the regularity 
of the solutions obtained is based on the use of the order completion 
method, first introduced and developed in [3]. 

As it happens, however, this order completion method reaches far be- 
yond the solution of systems of nonlinear PDEs, and in fact it can be 
applied to the solution of the surprisingly general equations 

(1.2) T(A) = F 
where 

(1.3) T : X — ► Y 
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is any mapping, X is any nonvoid set, while (Y, <) is a partially or- 
dered set, or in short, poset, while F e Y is given, and i 6 I is the 
sought after solution. 

Needless to say, in general, for a given F eY, there may not exist any 
solution A e X for the equation (1.2). Consequently, the setup (1.2), 
(1.3) may have to be extended. 

Customarily, such extensions assume suitable topologies on X and Y, 
and certain continuity properties for the mapping T in (1.3). 
However, as shown in [3], and also seen in the sequel, for the same 
purpose of solving the equations (1.2) in an extended setup, one can 
successfully use the order completion of the spaces X and Y. 

And one of the major advantages of such an approach is that such a 
method 

• does no longer differentiate between linear and nonlinear opera- 
tors T in (1.3), in case the spaces X and Y may happen to have 
a linear vector space structure, [3]. 

Two particularly convenient further features of the order completion 
method in solving general equations (1.2) are the following : 

• one obtains necessary and sufficient conditions for the existence 
of solutions, 

• one obtains explicit expressions of the solutions, whenever they 
exist. 

We shall present one of the general approaches resulting from the or- 
der completion method for solving equations of type (1.2). Further 
possible developments in this regard of the order completion method 
will be indicated. 



2. Pull-Back Order 

Without loss of generality, [3] , we shall assume that all the posets con- 
sidered are without a minimum or maximum element. Various notions 



3 



and results related to partial orders which are used in the sequel are 
presented in the Appendix. 

Given an equation (1.2), (1.3), we define on X the equivalence relation 

~t by 

(2.1) u^ T y T(u) = T(v) 

for u, v G X. In this way, by considering the quotient space 

(2.2) X T = X/^ T 

we obtain the injective mapping 

(2.3) T a : X T — ► Y 
defined by 

(2.4) X T 3U\ — ► T{u) e Y 

where u G U, that is, U is the equivalence class of u in Xt, while 
T(u) is defined by (1.3). 

At that stage, we can define a partial order <t on X T as being the 
pull-back by the mapping T~ in (2.3) of the given partial order < on 
Y, namely 

(2.5) U < T V T^{U) < T*(V) 

for U, V e Xt- The effect of the above construction is that we obtain 
the order isomorphic embedding, or in short OIE 

(2.6) T~ : X T — > Y 

As mentioned, without loss of generality we shall assume that the 
poset (A*, <t) has no minimum or maximum. 

And now, we consider the order completions X* and F # of X T , and 



respectively, Y. 

For simplicity, we shall denote by < the partial orders both on A* 
and Y#. In fact, as seen in (A. 20), these partial orders are the usual 
inclusion relations C among subsets of A, respectively, of Y. 

Then according to Proposition A.l in the Appendix, we obtain the 
commutative diagram of OIE-s 



Xn 



Y 



(2.7) 



X. 



c 



# 



c 



Y# 



Consequently, for U G X T and A e A*, we have in Y* the relations 

(2.8) T#(< U\) = <T^(U)} 

(2.9) T#(A) = (T^(A)) ul = sup y# { <T*(U)] \ U e A } 



3. Reformulation 

Now we can reformulate the problem of solving the general equations 
(1.2) as follows. Given F e Y#, find necessary and sufficient condi- 
tions for the existence of A e A*, such that 

(3.1) T*(A) = F 



4. Solution 



We note that (2.7) gives the inclusions 
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sup y# { T*(U) | U G X*, T#(U) C F } C 

C T#( sup x# {U\UeX*, T#(U) C F } ) C 

(4.1) 

C T#( inf x # {VlFel* F C T#(V) } ) C 
C inf y# { T # (V) \ V eX*, F C T # (\/) } 

Indeed, the first and last inclusions follow from Lemma A.l in the 
Appendix. As for the middle inclusion in (4.1), let U, V e X* be such 
that T#(C/) C F C T*(V). Then T#(C/) C T*(V), hence C/ < V, 
since T # is an OIE. It follows that 

sup x # { U | U e X*, T*(U) C F } < 

< inf x # { V | F G X*, F C T#(F) } 

and the proof of (4.1) is completed. 

We note further that the above inequality T*{U) C F C T # (V) also 
implies 

sup y# { T # (f/) | [/ G X#, T # (f/) C F } C F C 

(4.2) 

C inf y# { T#{V) \ V eX*, F C T # (\/) } 

Furthermore 

(4.3) { {/ e l| | T # (f/) C F } ^ 

since (A. 9) gives U — G X*, hence in view of (A. 4) - (A. 6) and 
(2.9) we have T#(U) = {T{(f)) ul = <p ul = ((f) u ) 1 = (Y*) 1 = C F. 

Returning now to the problem (3.1), we note that it is not trivial. 
Indeed, the OIE in (2.7), namely 

T* : X* — ► Y* 
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need not be surjective. Further T # need not preserve infima or suprema. 

However, in the next theorem we can obtain the following two general 
results : 

• a necessary and sufficient condition for the solvability of (9,12), 
and 

• the explicit expression of the solution, when it exists. 

Theorem 4.1. 

Given F eY*. 

1) The equation 
(4.4) T#(A) = F 

has a solution A G X*, if and only if, see (4.1) 

sup y# { T*{U) | U G X*, T*(U) C F } = 

(4.5) 

= inf y# { T#(V) | V e X*, F C T#(V) } 

2) This solution is unique, whenever it exists, see (2.7). 

3) When it exists, the unique solution A e X* is given by 

A = sup x# { U G X* | T # (f7) C F } = 

(4.6) 

= inf x # { V G X* | F C T#(y) } 

and, see (4.3) 

(4.7) { [/ G X# | T # (f/) C F }, { V G X* | F C T # (V) } ^ 
Proof 

From (4.4) follows that AeX*, T#(A) C F, thus 
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F = T#{A) C sup y# { T*(U) | U e X*, T#(U) C F } 
Similarly we have 

inf y# { T#(V) \ V eX*, F C T # (V) } C T # (A) = F 
thus (4.1) collapses to the seven equalities 

F = T#{A) = sup y# . . . = T#(sup x # . . . ) = 

= T#(inf x # . . . ) = inf y# . . . = T*(A) = F 
Thus in particular we obtain (4.5). 

The injectivity of T* will give (4.6), while (4.7) follows from (4.3) and 
the fact that we can take V = A. 

Conversely, let us assume (4.5). Then (4.1) collapses to the three 
equalities 

sup y # ... = T # (sup x # . . . ) = T # (inf x # . . . ) = inf y # . . . 

thus in view of the corresponding collapsed version of (4.2), we can 
extend the above three equalities to the following four 

sup y# ... = T # (sup x # . . . ) = T # (inf x # . . . ) = inf y # ... — F 

And now the injectivity of T # will give (4.4) and (4.6), while (4.7) 
follows as above. 

□ 

The above existence result is of a "local" nature, since it refers to a 
solution of the equation (3.1) for one given right hand term F E Y*. 
This result, however, can further be strengthened by the following 
"global" one which characterizes the solvability of (3.1) for all right 
hand terms F G Y* . Namely, we have, [3, pp. 190,191] 
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Theorem 4.2. 



The following are equivalent 



(4.8) 



T*(X*) D Y 



and 



(4.9) 



T#(X*) = Y* 



In each of these cases T # is an order isomorphism, or in short, 01, 
between X* and Y* . 

It is important to note the following two fact : 

• "Pull-back" type structures are customary when solving PDEs 
by functional analytic methods. Details in this regard are pre- 
sented in [3, chap. 12], while one well known classical example 
can be seen in section 7 in the sequel. 

• As shown in [3, chap. 13], one can consider in (2.7) far more 
general partial orders than the "pull-back" type ones, and still 
obtain solutions for nonlinear PDEs in (1.1) by the order com- 
pletion method. 



5. Applications to Nonlinear Systems of PDEs 

Let us now associate with a nonlinear PDE in (1.1) the corresponding 
nonlinear partial differential operator defined by the left hand side, 



(5.1) T{x,D)U{x) = F{x,U{x), . . . ,DPU{x), . . . ), x E tt 



Two facts about the nonlinear PDEs in (1.1) and the corresponding 
nonlinear partial differential operators T(x,D) in (5.1) are important 
and immediate : 



namely 
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• The operators T(x, D) can naturally be seen as acting in the 
classical context, namely 

(5.2) T(x, D) : C m (n) 3 U i — ► T(x, D)U G C°(Q) 

while, unfortunately on the other hand : 

• The mappings in this natural classical context (5.2) are typically 
not surjective even in the case of linear T(x,D), and they are 
even less so in the general nonlinear case of (1.1). 

In other words, linear or nonlinear PDEs in (1.1) typically cannot be 
expected to have classical solutions U G C m (Q), for arbitrary contin- 
uous right hand terms / G C°(Q), as illustrated by a variety of well 
known examples, some of them rather simple ones, see [3, chap. 6]. 
Furthermore, it can often happen that nonclassical solutions do have 
a major applicative interest, thus they have to be sought out beyond 
the confines of the classical framework in (5.2). 

This is, therefore, how we are led to the necessity to consider gen- 
eralized solutions U for PDEs like those in (1.1), that is, solutions 
U <£. C m (f2), which therefore are no longer classical. This means that 
the natural classical mappings (5.2) must in certain suitable ways be 
extended to commutative diagrams 

T(x,D) 

C m (tt) C°(fi) 



(5.3) 



c c 



X ► Y 

T 

with the generalized solutions now being found as 
(5.4) U eX\ C m (Q) 

instead of the classical ones U G C m (fl) which may easily fail to exist. 
A further important point is that one expects to reestablish certain 
kind of surjectivity type properties typically missing in (5.2), at least 
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such as for instance 



(5.5) C°(Q) C T(X) 

As it turns out, when constructing extensions of (5.2) given by com- 
mutative diagrams (5.3), we shall be interested in the following some- 
what larger spaces of piecewise smooth functions. For any integer 
< / < oo, we define 

/r r\ ni /n\ f o to 3 T C closed, nowhere dense : 
(5-6) CM = [u-M^R ueC ^ T) 

and as an immediate strengthening of (5.2), we obviously obtain 



(5.7) T(x, D) C™ d (Sl) C CM 

The solution of the nonlinear PDEs in (1.1) through the order com- 
pletion method will come from the construction of specific instances 
of the commutative diagrams (5.3), given by 



(5.8) M%(£1) 



T(x,D) 



T 



bijective 



where the operation ( ) # means the order completion, [3], of the re- 
spective spaces, as well as the extension to such order completions of 
the respective mappings, see (2.7). It follows that in terms of (5.3), 
we have 



x = M%(n)*, Y = M (n)*, T = T* 
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thus we shall obtain for the nonlinear PDEs in (1.1) generalized solu- 
tions 

(5.9) U e M^(n)# 

Furthermore, instead of the surjectivity condition (5.5), we shall at 
least have the following stronger one 

(5.10) CM C T*(M?(m 

So far about the main ideas related to the existence of solutions of gen- 
eral nonlinear PDEs of the form (1.1). Further details can be found 
in [3,1,4-6]. 

As for the regularity of such solutions, we recall that, as shown in [1], 
one has the inclusions 

(5.11) M°(Vl)# C MesiSl) 

where Mes (Q) denotes the set of Lebesgue measurable functions on f2. 
In this way, in view of (5.8) and (5.9), one can assimilate the general- 
ized solutions U of the nonlinear PDEs in (1.1) with usual measurable 
functions in Mes (Q). 

Recently, however, based on results in [1,4-6], it was shown that in- 
stead of (5.11), one has the much stronger regularity property 

(5.12) M°{Q)# C M{tt) 

where H (Q) denotes the set of Hausdorff continuous functions on Q. 
Consequently, now one can significantly improve on the earlier regu- 
larity result, as one can assimilate the generalized solutions U of the 
nonlinear PDEs in (1.1) with usual functions in H(fi). 

Regarding systems of nonlinear PDEs such as in (1.1), with possibly 
associated initial and/or boundary value problems, it was shown in [3, 
chap. 8] the way they can be dealt with the above order completion 
method. 
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In this respect, a surprising advantage of the order completion method 
is the ease, when compared with the usual functional analytic ap- 
proaches, in dealing with initial and/or boundary value problems. 

6. Beyond "Pull-Back" Partial Orders 

As mentioned at the end of section 4, and presented in full detail in 
[3, chap. 13], the order completion method in solving large classes 
of nonlinear systems of PDEs of the type (1.1) is not limited to the 
use of "pull-back" type partial orders in (2.7), (5.3) and (5.8). In 
fact, a large class of more general partial orders can be defined on 
the domains «M™(fi) of the respective PDEs, and stil obtain for them 
solutions in the corresponding order completions. 

7. Use of "Pull-Back" in Functional Analytic Solution Meth- 
ods 

As presented in detail in [3, chap. 12], functional analytic methods 
used for solving PDEs do often employ topologies obtained by "pull- 
back" . Here we present shortly one of the classical such examples. Let 
us consider on a bounded Euclidean domain Q, which has a smooth 
boundary d Q, the following familiar linear boundary value problem, 
usually called the Poisson Problem 

A U(x) = f(x), xen 

(7.1) 

U = on dQ 

As is well known, for every given / e C°°(Q), where fl denotes the 
closure of Q, this problem has a unique solution U in the space 

(7.2) X = { v G C°°(n) | v = on d Q } 
It follows that the mapping 

(7.3) X3v i — > || Ad || x,2 ( n ) 
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defines a norm on the vector space X. Now let 

(7.4) Y = 

be endowed with the topology induced by L 2 (Q). Then in view of 
(7.1) - (7.4), the mapping 

(7.5) A : X -> Y 

is a uniform continuous linear bijection. Therefore, it can be extended 
in a unique manner to an isomorphism of Banach spaces 

(7.6) A : X -> Y = L 2 {Q) 

In this way one has the classical existence and uniqueness result 
V feL\Q) : 

(7.7) 3! [/el : 

AU = f 

The power and simplicity - based on linearity and topological comple- 
tion of uniform spaces - of the above classical existence and unique- 
ness result is obvious. This power is illustrated by the fact that the 
set Y = L 2 (Q) in which the right hand terms / in (7.1) can now be 
chosen is much larger than the original Y = C°°(Q). Furthermore, 
the existence and uniqueness result in (7.7) does not need the a priori 
knowledge of the structure of the elements U G X, that is, of the 
respective generalized solutions. This structure which gives the regu- 
larity properties of such solutions can be obtained by a further detailed 
study of the respective differential operators defining the PDEs under 
consideration, in this case, the Laplacean A. And in the above specific 
instance we obtain 

(7.8) X = H 2 (Q)nH^(Q) 

As seen above, typically for the functional analytic methods, the gen- 
eralized solutions are obtained in topological completions of vector 
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spaces of usual functions. And such completions, like for instance the 
various Sobolev spaces, are defined by certain linear partial differential 
operators which may happen to depend on the PDEs under consider- 
ation. 

In the above example, for instance, the topology on the space X ob- 
viously depends on the specific PDE in (7.1). Thus the topological 
completion X in which the generalized solutions U are found accord- 
ing to (7.7), does again depend on the respective PDE. 



Appendix 

We shortly present several notions and results used above. A related 
full presentation can be found in [3, Appendix, pp. 391-420]. 

Let (X, <) be a nonvoid poset without minimum or maximum. For 
a G X we denote 

(A.l) < a] = {x E X | x < a}, [a >= {x E X | x > a} 

We define the mappings 

(A.2) X D A i > A u = f] aeA [a > C X 

(A.3) X D A i > A 1 — f] aeA < a] C X 

then for A C X we have 

(A.4) A u = X A 1 = X A = (j) 

(A. 5) A u = <fi A unbounded from above 

(A. 6) A 1 = (f) <^=^ A unbounded from below 

Definition A.l. 

We call A C X a cut, if and only if 
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(A.7) A ul = A 
and denote 

(A.8) X* = {A C X | A is a cut} C P(X) 

□ 

Clearly, (A.4) - (A. 6) imply 

(A.9) (f), x ex* 

therefore 

(A.10) X#^0 

Given A, B C X, we have 

(A.ll) A C £ =>■ A M D B u , A 1 D B l 

(A.12) A C A ul , A C 

(A. 13) A"'" = A", A M = A' 

Consequently 

V ACI : 

*) A ui G X # 

(A.14) * *) V B E X* : 

AC5 =>■ A ul C £ 

B C A =>• B C A"' 

therefore 
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(A.15) X* = {A ul \ACX} 
Given x G X , we have 

(A. 16) {x} u = [x >, {x} 1 =< x], [x > l =< x], < x] u = [x > 

(A.17) {x} ul =<x\, {x} lu =[x> 

We denote for short 

= x « j = x i^ { x ya = x ui^ { x y u = x iu^ _ _ _ 

Given A G X # , we have 

/ 3 : \ 

(A.18) <^A^X^ 

\ <a] C A C <b] J 
We shall use the embedding 

(A.19) X 9 x ^ x u/ = x' =< x] G X* 

We define on X* the partial order 

(A.20) A < B ACB 

Definition 2.1. 

Given two posets (X, <), (Y, <) and a mapping </? : X — > Y. We call 
ip an order isomorphic embedding, or in short, OIE, if and only if it is 
injective, and furthermore, for a, b G X we have 

a < b ip(a) < (p(b) 

An OIE <£> is an order isomorphism, or in short, 01, if and only if it is 
bijective. 

□ 
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The main result concerning order completion is given in, [2] : 

Theorem ( H M MacNeille, 1937 ) 

1) The poset <) is order complete. 

2) The embedding X — X* in (A. 19) preserves infima and suprema, 
and it is an order isomorphic embedding, or OIE. 

3) For A E X*, we have the order density property of X in X*, 
namely 

A = sup x # {x l | x E X, x l C A} = 

(A.21) 

= inf x # {x l | x G X, AC x 1 } 

□ 

For A C X, we have 

(A.22) A ul = sup x# {x l | x e A} 

Given Ai e X # , with i e we have with the partial order in X # the 
relations 

(A.23) sup JG/ Ai = inf {A e X# \ {J^AQA} = ([j^A^ 1 

M ieI A, = sup {A eX* | A C f| ie/ A} = (fW^ = 

(A.24) 

= rw^ 

Extending mappings to order completions 

Let (X, <), (Y, <) be two posets without minimum or maximum, and 
let 

(A.25) ip:X — > Y 
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be any mapping. Our interest is to obtain an extension 

v # : x * — > r# 

For that, we first extend </? to a larger domain, as follows 
(A.26) y? # ■ V(X) — ► Y* 
where for ACIwe define 

(A.27) V #(A) = (ip(A)) ul = sup Y# {< tp(x) ] \ x <E A} 

and for any mapping in (A. 25), we obtain the commutative diagram 



X 3 x \- 



<f(x) e Y 



(A.28) 



V(X) 3 {x}\- 



# 



ip#( x ) = < ip{x) \ e Y* 



Proposition A.l. 

1) The mapping : V(X) — >■ in (A. 36) is increasing, if on 
V(X) we take the partial order defined by the usual inclusion "C". 

2) If the mapping tp : X — > Y in (A. 35) is increasing, then the 
mapping tp# : V(X) — >■ in (A. 36) is an extension of it to X*, 
namely, we have the commutative diagram 



X 3 x h 



<f(x) e Y 



(A29) 

X* 3 < x] h 



# 



v? #(< x ] ) = < ip(x)} e Y* 
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3) If the mapping (p : X — ► Y in (A. 25) is an OIE, then the mapping 
(p# : V(X) — > y* in (A. 26) when restricted to X#, that is 

(A.30) if* : X* — ► Y# 

as in (A. 29), is also an OIE. 



Lemma A.l. 

Let in general // : M — > N be an increasing mapping between two 
order complete posets, then for nonvoid £CMwe have 

(A. 31) ^(inffcf E) < mix fi(E) < sup N n{E) < /i(sup M E) 

Proof 

Indeed, let a = inf^ E G M. Then a < b, with b G E. Hence 
fi{a) < with b <E E. Thus /x(a) < infjv ^(E), and the first 

inequality is proved. 

The last inequality is obtained in a similar manner, while the middle 
inequality is trivial. 

□ 



References 

[1] Anguelov R, Rosinger E E : Hausdorff continuous solutions 
of nonlinear PDEs through the order completion method. 
Quaestiones Mathematicae, Vol. 28, 2005, 1-15, arXiv : 
|math.AP/0406517| 

[2] Luxemburg W A J, Zaanen A C : Riesz Spaces, I. North-Holland, 
Amsterdam, 1971 



20 



[3] Oberguggenberger M B, Rosinger E E : Solution of Continuous 
Nonlinear PDEs through Order Completion. Mathematics Stud- 
ies VOL 181, North-Holland, Amsterdam, 1994 

[4] Rosinger E E : Hausdorff continuous solutions of arbitrary con- 
tinuous nonlinear PDEs through the order completion method. 
|math.AP/0405546| 

[5] Rosinger E E : Can there be a general nonlinear PDE theory for 
the existence of solutions ? ||math.AP/ 0407026 

[6] Rosinger E E : Solving large classes of nonlinear systems of PDEs. 
|math.AP/0505674| 



21 



